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Abstract 

^ , We seek to complement Nelson's work on the two-slit experiment 

by showing that the two-slit process, whose density exhibits the char- 
acteristic interference pattern, may be obtained as the model after 
^ , the beam has reached the screen by means of a variational mecha- 

?—i ' nism. The one-slit process, modeling the beam before it reaches the 

screen, plays the role of a reference model. 
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1 Introduction 

When a beam of electrons originating from a source S goes through a crystal, 
one observes a diffraction pattern resembling that characteristic of interfer- 
ence in wave motion. In a celebrated gedanken experiment [1, pp. 2-5], the 
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crystal is replaced by a screen with two slits. As it is well-known, according 
to quantum mechanics, this is an instance of a superposition of two states. 
The two possibilities of going through the top or the bottom slit interfere. 
Hence, the probability of arrival at the second screen is not simply the sum 
of the probabilities of coming through the top and the bottom slit. 

In [2, 3], Edward Nelson made a remarkable calculation showing that, 
according to his version of stochastic mechanics, the observed interference 
effect can be explained within the frame of classical probability. In particular, 
the probability of arrival is indeed the sum of the probabilities of coming 
through the top and the bottom slit. 

In this note we seek to complement what Nelson did. We model the emis- 
sion from the source through a suitable Gaussian process (one-slit process). 
We then take the latter as a reference process in a variational problem that 
takes into account the presence of the screen with the two slits. The solution 
of the variational problem, suitably extending on the results of [4, 5], is the 
stochastic process that Nelson has analyzed (two-slit process) whose density 
profile is familiar in wave interference. 

2 Elements of the Nelson-Guerra stochastic 
mechanics 

Nelson's stochastic mechanics is a quantization procedure for classical dy- 
namical systems based on diffusion processes. Following some early work 
by Feynes [6] and others, Nelson and Guerra elaborated a clean formula- 
tion starting from 1966 [7, 8, 9] and references therein, showing that the 
Schrodinger equation could be derived from a continuity type equation plus a 
Newton type law, provided one accepted a certain definition for the stochastic 
acceleration. The Newton-Nelson law was later shown to follow, in analogy 
to classical mechanics, from a Hamilton-like stochastic variational principle 
[10, 11]. Other versions of the variational principle have been proposed in 
[2,12,13,14]. 

Consider the case of a nonrelativistic particle of mass m. Let {ip(x, t); t < 
t < satisfying the Schrodinger equation 
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be such that 



IIW>|||e Lj, c [t ,+oo) 



(2) 



This is Carlen's finite action condition. Under these hypotheses, the Nelson 
measure P may be constructed on path space, [15], [16], [12, Chapter IV], and 
references therein. Namely, letting f2 := C([t , ti], R n ) the n-dimensional con- 
tinuous functions on [io,£i], under the probability measure P, the canonical 
coordinate process x(t,u) = oo(t) is an n-dimensional, Markov, finite-energy 
diffusion process {x(t); t < t < t{\ , called Nelson's process, having (forward) 
Ito differential 



dx(t) 



h 



m 



V (3« logi>(x(t), t) + 3 log if>(x(t), t)) 



dt + \ —dw(t), (3) 
m 



where w is a standard, n-dimensional Wiener process. Moreover, the proba- 
bility density p(-,t) of x(t) satisfies Born's relation 



p(x,t) 



(x,t)| 2 , Vfe[f ,*i]. 



(4) 



3 Nelson's treatment of the two-slit experi- 
ment 

We recall some essential points in Nelson's analysis in [2, 3]. Let x G M 3 , and 
consider free motions, i.e. V = 0, in a frame of reference comoving with the 
beam. For the sake of simplicity, take m — 1, h — 1 and to — 0. The slits 
are located at ±a, where a is a vector in R 3 . Let 

a\ 1/4 , \x\\ 



= - exp 



nj " v 2A ' 

be the initial condition at time to- Then, the corresponding solution of 

dip a , / _ 

if = 2^'- (5) 



IS 

Wl , t)= gy / \ A+it) -v 2exp /_M_v (6) 
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Nelson calls the one-slit process the Gaussian process associated to this evo- 
lution whose forward and backward drifts are given by 

He then takes as initial condition for (1) the vector 

V'i(x) = 7(^0(2; - a) + i>o(x + a)) , (7) 

where 7 is a constant (close to l/y/2 if \a\ is large with respect to A) to make 
■01 a unit vector in L 2 . By the linearity of the Schrodinger equation, the 
corresponding solution is 

ipi(x, t) = 7 (V>o(z - a, t) + ^o(z + a, t)) . (8) 

Nelson goes ahead and computes the forward drift of the two-slit process 
corresponding to As only the direction joining the slits is of interest, 
the problem can be reduced to a one-dimensional problem (hence, from now 
on, a,x £ K). He finds that, for small times, the drift is nearly the same 
as for the one-slit process. For larger times, however, the drift of the two- 
slit process becomes enormous in certain regions, repelling particles from 
there. As a result, the probability density of the two-slit process after some 
time exhibits alternate peaks and valleys resembling those observed in wave 
interference. 



4 A stochastic variational principle 

Let v be the velocity of the beam, let d be the distance between the source 
and the screen with the two slits. Let T = d/v be the time it takes the beam 
to reach the screen. In order to keep t — the time when the beam is at the 
screen, we assume that the beam is emitted from the source S at time — T. 
The one-slit process with A very small associated to ip (x,T + t) as in (6) 
models effectively the beam of particles between times — T and 0. Suppose 
now we impose at time t — the probability density 

Po( x ) = 2 {\^o(x - a)\ 2 + \ip (x + a) | 2 ) . 

We define Xq to be the family of real- valued stochastic processes {x(t); —T < 
t < 0} with continuous paths and satisfying the following properties: 
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1. x(t) has a nowhere vanishing probability density p(x,t) of class C 1 for 
all t G [— T, 0] with marginal probability density po a t time 0; 

2. £{/° T |J>gp(x(t),t)| 2 dt}<oc; 

3. there exists for each t G [— T, 0] a random variable v(t), called the 
current drift, such that, if f(x,t) is of class C 2 ' 1 satisfying 

< oo, 



E{J_ T \-^-f(x(t),t)\ 2 dt} 



we have 



Notice that property (iii) is satisfied by a rather large class of processes 
including differentiable processes, finite-energy processes with constant dif- 
fusion coefficient [17, 2]. It also includes Markovian diffusion processes with 
local diffusion coefficient given the results in [21, 20, 19, 18]. The current 
drift v is just the semi-sum of the forward and backward drifts. Let V denote 
the family of finite-energy stochastic processes on [— T, 0] . In the following 
variational problem, we take one-slit process {x r (t + T); — T < t < 0} with 
small A associated to {ipo(x,T + 1); —T < t < 0} as reference model on the 
time interval [— T, 0]. Write ipo(x, T + t) = p r (x, t) exp[^S r (x, t)]. Rather 
than setting up a variational problem with complex-valued drifts of the form 
v — iu as in Section VIII in [4] and [5], we formulate here a problem with 
distinct real- valued velocities v and u as in Section IV in [4] and [22]. For 
(x,v,u f ) G (X , V, V), define the functional 



I(x,v,u) -=E\f_ 



-0 
T 



-£( t,,(t) -s iogftWt) ' i)> 



dt- S r (x(0),0) 



Notice that the integrand is the same as in the Guerra-Morato action [11] 
where differences of drifts replace drifts. We consider the stochastic differen- 
tial game 
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min max min I(x,v, u') 
x&x vav u'ev 

subject to 

v(-) is the current drift of x, 
u'(t) = log p(x(t), t) Vt , 

where p(-,t) is the probability density of x(t). We say that (x*,v*,u'*) 
is a saddle-point equilibrium solution of the game if for all (xi,v*,u') and 
(x2,v,u'*) in Xq x V x V satisfying the constraints we have 

I(x!,v*,u') < I(x*,v*,u'*) < I(x 2 ,v,u'*). 

Let F : (M x [io,*i]) — > R be of class C 1 satisfying the condition in property 
hi). Let A : (R x [t Q , ti\) — > R be of class C 1 . For such a pair (F, A), define 



A F ' A (x, «') = F {F(a;(0), 0) - F(x(-T), -T) 



d 



u'(t)-—\ogp(x(t),t) 



dt 



Observe that for all triples (x, v,u') in X x V x V satisfying the constraints, 
we have A F ' A (:r, t>, -u') = 0. Obviously, if (x,v,u') E (X x V x V) satisfying 
the constraints is a saddle-point solution for (7 + A F,A ), then it also solves 
the original problem with cost function /. 

It is possible to rewrite A F ' A in a form more suited for our purposes using 
argument similar to that in [11]. Observe that, for any x G X , we have 

E^X(x(t),t)^logp(x(t),t)^ = J Rn \(x,t)-^p(x,t)dx 



d_ 

1" dx 



\(x, t)p(x, t)dx 



- E {-?-\(x(t),t)\, 



[ dx 



where, in the integration by parts, we have used the natural boundary con- 
dition at infinity for p(-,t). Thus, our A F,A functionals now have the form 



A F ' x (x, v, u') = E{F(x(0), 0) - F(x(-T), -T) 



dF , \ d 



d 



{x(t),t) - X(x(t),t) ■ u\t) - —\(x{t),t) dt}. 



The variational analysis now follows the same lines as in [22]. Consider now 
the unconstrained problem 

min minmaxfJ + A F,x )(x, v, u'). 
xex vev u'ev 

For each fixed x G X , and each t G [— T, 0], we study the finite- dimensional 
problem 



minmax < -m(v — S r (x(t), t)) 2 (V — — — log p r (x(t), t)) 2 

«eRu'eR\2 V mdx rK Kh JJ 8m K dx BH K w ' ;; 

OF d d 1 

- — (x(t),t) - v—F(x(t),t) + X(x(t),t)u' + —\(x(t),t)\ . 

We get the optimality conditions 

4:171 

U'* x {t) = —\0gp r (x(t),t)) + -j^\(x(t),t). 

Remark 4.1 If a stochastic process with the prescribed v*(t) and u'*(t) does 
exist, then the first optimality condition implies that it is a Markov process 
and the second that A is given by 

^i 1 * (£)<*■*>• 

where p* x (x,t) denotes the probability density of x*(t). 

Notice that v x and u' x belong to V. Consider next the minimization of 

(/ + A F ' A )(x, <*) = E {F(x(0), 0) - F(x(-T), -T) 



L 



+ 

-T 



*2m d d 

-^\(x(t),t) 2 + —\og Pr (x(t),t)) ■ \(x(t),t) + — ■ \(x(t),t) 
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on the space Xq. We wish to choose F and A such that the functional becomes 
constant with respect to the process x e X . Suppose that the pair (F, A) 
satisfies onlx [— T, 0] the equation 

— —(—fX —-S — F-— A 2 - — lo A-— A-0 (9) 

dt 2m \ dx J mdx r dx h 2 dx ^ ^' dx 

and the boundary condition F(x, t ) = —S r (x, 0). With this choice of (F, A), 
(I + A F > x )(x,v*,u'*)) = £'{F(x(0), 0)} which is constant on X (all processes 
have the same marginal density at t — 0). Hence, any x in X solves with 
(v x ,u' x ), the unconstrained problem. 



5 Solution to the variational problem 

In view of Remark 4.1, we now write 



for some scalar C 1 function G. Equation (9) then becomes 

d d 

0. 



dF 1 / d \ 2 1 d „ d „ h 2 



dt ^ 2m \ dx^j ~^mdx^ r dx^ 2m 



'<L G X + —G + — lo —G 
dx i dx 2 dx dx 



(10) 

In the ideal case where -j^S r = and p r = const, this equation reduces to one 
of the Madelung equations. If F and G may be found satisfying (10) with 
F(x, to) = —S r (x, 0), then any x in X Px solves the unconstrained problem. If 
we can find one x* in X pi that also satisfies the constraints 

Id Id 
v*(t) = — —S r (x(t),t)) H —F(x(t),t) is the current drift of x, 

mdx mdx 

d d d 

— log p*(x(t),t) = —logp r (x(t),t) + 2—G(x,t)Vt, 

then it solves the original stochastic differential game. In that case, we define 

S*(x,t) := S r (x,t) + F(x,t), (11) 
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so that v*(t) = -^-^S*(x(t),t)) and we have 



G(x,t) = ±]og&(x,t). (12) 
z p r 



Suppose that we have found such a process x* with density p*. As observed 
before, it is a Markov process and so is the Nelson reference process x r . Then 
the corresponding Fokker-Planck (continuity) equations read 

dp* | did = Q 
<9t dx mdx ' 

£ + < l4 » 

Using these equations, and relations (11)-(12), we get that F and G must 
satisfy another equation 



dG _ <9(|logf) _ 1 



1 13,19. 

— tt > + ~ir(-ir s rPr 

p* ox m ox p r ox m ox 



dt dt 2 



;(S* - S r ) -^—S*—— bg ^ — (5* - S r )—\og Pr 



2m dx 2 2m dx dx p r 2m dx dx 

——F--—F—G+-—S—G-——F—\o 

2m dx 2 mdx dx mdx r dx 2m dx dx ^ ^ r 

Thus, F and G also satisfy 

— + -— F— G+— — F+— — F — lo - —— S— G-Q (15) 

dt m dx dx 2m dx 2 2m dx dx mdx r dx 

In the ideal case where -§^S r = and p r = const, this equation reduces to 
the other Madelung equation. Moreover, 

G(s,0) = £log^(s,0). (16) 
z p r 

Define 9 = exp(G + j-F). Then, equations (10)-(15) imply that 9 satisfies 

90_ + (]_d_ Sr + A A* _ *L&_ e = (17) 

<9t \mdx r im dx ') dx 2m dx 2 

Define now ip n : = (p*) 1/2 exp jrS* = ip r 9. It follows from (5) and (17) that 
ip n also satisfies (5) and ip n (x,0) = (po(x)) 1 ^ 2 . 
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6 Conclusion 



Hence, we find that the solution to the stochastic differential game is the 
Nelson process associated to another solution of the same Schrodinger equa- 
tion (5), with the probability density po at time and S*(x, 0) = (so that 
v(0) = 0). In other words, the variational principle shows that the process 
with density p at time and v(0) = that replaces the one-slit reference 
process is precisely the two-slit process studied by Nelson in [2, 3]. 

As observed by Nelson in [3], the probability density po{x) can, for all 
practical purposes, be replaced by \ip\{x, 0)| 2 , where ipi was defined in (8). 
Thus, the new process, obtained from the variational principle, will exhibit, 
after a time sufficiently large with respect to a, the typical interference pat- 
tern in the probability density. 

We have made here no attempt at a physical interpretation of the varia- 
tional mechanism. We refer the reader to [23, 24, 25], and references therein, 
for a discussion of various crucial topics in stochastic mechanics from the 
physical standpoint. 
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